Thermally generated vortices, gauge invariance and electron spectral function in the 

pseudo-gap regime 
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Starting from classical vortex fluctuation picture, we study the single electron properties in the 
pseudogap regime. We show that it is the gauge invariant Green function of spinon which is directly 
related to ARPES data in the pseudogap regime instead of the non-gauge invariant one. We find 
that the random gauge field from the thermally generated vortices completely destroys the coherent 
spinon motion and leads to collective excitations pertinent to non-Fermi liquid behaviors. The 
Energy Distribution Curves (EDC) show broad peaks, while the Momentum Distribution Curve 
(MDC) show rather sharp peaks with Lorenz form. The local density of state at zero energy scales 
as the inverse of Kosterlize-Thouless length. These results are qualitatively consistent with the 
ARPES data in the pseudo-gap regime. 



Although the superconducting state of high T c 
cuprates can be described approximately as pair con- 
densates with a d-wave symmetry and well-defined 
quasi-particle excitations above them, the anomalous 
properties of the normal state are still poorly under- 
stood. Many experiments such as Angle-Resolved Pho- 
toemission (ARPES) M, Scanning Tunneling Microscopy 
(STM) 0, Inelastic Neutron Scattering experiments g] 
and transport characteristics provided evidence for 
the opening of a gap in the electronic excitation spectrum 
above the critical temperature T c up to some characteris- 
tic temperature T* in the underdoped cuprates, this en- 
ergy gap is dubbed " pseudo-gap" ( Fig.l). Both ARPES 
HJ and STM have indicated that in the pseudo-gap 
regime, the pseudogap possesses the same d x 2_ y 2 sym- 
metry as the superconducting gap below T c . ARPES de- 
tected sharp quasi-particle peak below T c , but only broad 
spectral weight above T c . 
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Fig 1: There is a quantum phase transition from d-wave super- 
conductor to the underdoped regime at x = x c . The Quantum Crit- 
ical (QC), Quantum Disordered (QD) and Vortex plasma regimes 
are denned in Jl4| . There is no spin-charge separation in QD and 
QC regimes. But there is approximate spin-charge separation in 
the vortex plasma regime (see the text). 

In this paper, we are trying to explain the ARPES data 



near the four nodes (±0.397r, ±0.397r) M from the phys- 
ical picture proposed by Emery and Kivelson and others 
1^-^: the finite temperature superconducting transition 
is a Kosterlize-Thouless (KT) transition driven by vor- 
tex pair un-binding, namely T c — Tkt- Below Tkt, 
the phase fluctuation is small, the positive and nega- 
tive vortex pairs are tightly bound together, above Tkt, 
the long-range superconducting coherence is destroyed 
by strong phase fluctuation which leads to the gener- 
ation of free vortices with both parity, but the ampli- 
tude of the local superconducting order parameter re- 
mains non- vanishing up to T* . The quasi-particles in 
the pseudo-gap regime are the remnants of those in the 
superconducting phase and are strongly scattered by the 
free vortices above Tkt- We are trying to calculate the 
real electron spectral weight due to these scatterings and 
compare with the experiments. 

The quasi-particle scattering due to the Volovik effect 
fl0| in the pseudogap re gim e was discussed in Refs. 0. 
As first stressed in Ref. Jl^] , there is also Aharonov and 
Bohm (AB) phase effect from hc/2e vortices. Although 
the AB phase effect only provides a periodic potential 
in the vortex lattice states discussed in |IJ,[r|], it starts 
to scatter quasi-particles when the vortices are randomly 
distributed. As pointed out in [^5|, the elastic scattering 
due to the AB phase from the disordered vortex array 
dominate over that from the Volovik effect. In Mil- 
lis and the author studied the zero temperature quantum 
phase transition from d-wave superconductor to the un- 
derdoped regime driven by the condensation of quantum 
generated hc/2e vortices. In this paper, we apply the 
method developed in Refs. to study the electron 

spectral weight in the classical vortex fluctuation regime 
( the vortex plasma regime in Fig.l ). We derive the im- 
portant Eq.^ by two approaches. First, replacing vortex 
density n v generated by external magnetic field in Jig] by 
the free vortex density n/(T) above T > Tkt, we lead to 
Eq.^| directly from the results in Ref. |lq] . Second, start- 
ing from the zero temperature effective action developed 
in Ref. jl6| , we show that the action reduces to EqJ^ after 
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treating the vortices classically. We show that the ran- 
dom gauge field a a in Eq.^ from the thermally generated 
vortices completely destroys the coherent quasi-particle 
picture and leads to collective excitations pertinent to 
Non-Fermi liquid character. The Energy Distribution 
Curves (EDC) ( at fixed k ) show broad peaks, while 
the Momentum Distribution Curve (MDC) ( at fixed ui ) 
show rather sharp peaks with Lorenz form. Volovik effect 
leads to finite local DOS scaled as the inverse of the KT 
correlation length. These results are qualitatively con- 
sistent with the ARPES data in the pseudo-gap regime. 
The ARPES spectrum in QC and QD regimes will be 



discussed elsewhere 24 



In Ref. [[L5|, starting from BCS Hamiltonian, perform- 
ing the singular gauge transformation and expanding 
H s around the node 1 where p = (p^,0), Ye obtained 
H s = Hi + H c where the linearized Hamiltonian Hi is 
given by 0: 

Hi = v f ( Px + at)r z + v A (p y + a^T 1 + v f v x (r) 
+ (1^2,x^y) (1) 

where v s = \{vf+vf) = |V0—|A is the total superfluid 
momentum and = \{v^ — — \ (y<f>A — V^>b) is the 
internal gauge field. We get the corresponding expression 
at node 1 and 2 by changing Vf —> —Vf,v<\ — > — va in the 
above Eq. For explicit spin SU (2) invariant formulation, 
see Ref. @. 

The curvature term H c can be written as: 



H c = — [{n Q ,t> Q } 
m 



n 2 



:r 3 + |i{n II njr 1 ] (2) 



Where ft = p + is the covariant derivative. H c takes 
the same form for all the four nodes |ll| . 

Just like at filling factors at v — 1/3 or v = 1/2 QH 
system |2(J, the transport properties can be studied di- 
rectly in the transformed Hamiltonian. But the single 
particle electron Green function is much more difficult to 
calculate [ pl[ . In the present problem, the general form 
of electron annihilation operator was given in Ref. jl6| : 

C a (x)= ]T[e^V^V/ Q ^ia 



e a pe 



(3) 



where i is the node index, 1, 2 are p — h indices and a, /3 
are spin indices. 

All the above equations can be applied straightfor- 
wardly to the pseudogap regime above T > Tkt where 
vortices are unbound. There are equal number of free 
positive and negative vortices generated by thermal fluc- 
tuation. On the average, the vanishing of and are 
automatically ensured. Furthermore, the thermal gen- 
erated vortices are moving around and their smearing 
effect makes the U(V) character stronger than their Zi 



character in randomly-pinned static vortex array case. 
Replacing vortex density n v generated by external mag- 
netic field in pi by the free vortex density n / (T) above 
T > Tkt, we can describe the effective action by spinons 
moving in a static magnetic field ( there is no feedback 
effects on the gauge field propagator from the fermions ) 
ppl: 
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£ = i>b»(d fl -iai)i> a + —(f% ) 



(4) 



Where a, = 1,2 are the space components only; 
a = 1, • < - , JV = 4 are 4 species of Dirac fermion and 
n f {T)~ ^ 2 {T) where £(T) is the KT correlation length 
at T > T c . 

In Ref. the authors derived an effective action 

describing the zero temperature quantum phase transi- 
tion from d-wave superconductor to underdoped regime 
around the critical doping x c . At zero temperature, the 
vortices have to be treated quantum mechanically, one 
has to perform two singular transformations which are 
dual to each other to quasi-particlcs and moving vor- 
tices respectively to keep all the possible commutation 
relations intact. Just like conventional singular gauge 
transformation pp| leads to conventional Chern-Simon 
(CS) term, the two mutual singular gauge transforma- 
tions lead to a mutual CS term. In the effective action, 
quantum generated vortices couple to quasi-particles by 
the mutual CS term, the vortices are also interacting with 
a charge U(l) gauge field a M mediating the long-range 
current- current logarithmic interactions between the vor- 
tices. The U(l) charge fluctuation leads to mass terms 
not only for v, but also for the spatial component a^. 
There is no gapless dynamic gauge fluctuations in the 
Cooper-pair picture. The precise nature of the zero tem- 
perature quantum critical point (QCP) between the d- 
wave superconductor and some unknown charge ordered 
state is still under investigation. The electron spectral 
functions in the QC and QD regimes are controlled by 
this QCP and will be presented elsewhere p4| . 

However, the spectral function in the vortex plasma 
regime around Tkt in Fig.l are determined by the clas- 
sical hydro-dynamics of vortices which is the focus of the 
present paper. The vortices being treated classically [p^[ , 
their commutation relations can be neglected, the dual 
singular gauge transformation in the vortices is not nec- 
essary, the Berry phase term for the boson ( the linear 
derivative term) can also be neglected [jl9]]. More specif- 
ically, only the time component of the charge gauge field 
ao is kept to mediate the long-range density-density log- 
arithmic interaction between the vortices. Obviously, ao 
couples to the vortices the same way as the supcrfluid 
velocity v a couples to the spinon, its fluctuation leads to 
a mass term for the v — v correlation. However, being 
orthogonal to the spatial component a 



a v corre- 



lator remains gapless and is given by the Maxwell term 
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in Eq.[|. This is exactly consistent with the results first 
derived in 

It is important to point out that what measured by 
ARPES is the spectral function of the real electron in 
Eq.^ which carry both spin and charge, instead of that 
of the spinon ip which carry spin only. It is vital to 
emphasize that the electron Green function G(x,t) =< 
C a (0, 0)C|(s, t) > is gauge invariant under the internal 
gauge transformation but that of the spinon ip is 
not ! We make the following gauge invariant mean field 
approximation for the electron Green function: 



G(x,t) =< e ^)/ 2 e-^°)/ 2 > 



+ e 



< ^ la (x,t)e i r^^ d ^L(0,0) > 
' f <4(^)e-I'"* (l) V(0,0) 



> 



(5) 



As shown in Ref. Jl5|, it is vital to make a mean field 
approximation to respect the exact T symmetry. Very 
similarly, in the pseudo-gap regime, it is equally impor- 
tant to make a mean field approximation to respect the 
exact gauge invariance. Any approximation violating the 
gauge invariance will lead to completely wrong conclu- 
sions ! The vortices being treated classically, only the 
lj = component is kept. So the vortices fluctuate much 
slowly than the fermions which are treated completely 
quantum mechanically. We expect the the spin-charge 
separation indicated in EqJ^ is a good approximation in 
the vortex hydro-dynamic regime. However, as shown in 
p6| , there is no spin-charge separation at T = 0, in QC 
and QD regimes (Fig.l). This is in contrast to Z-i gauge 
theory j |25| ] where there is true spin-charge separation at 
T = 0. 

The first factor in Eq.|| is the classical correlation func- 
tion of the vortices, at T > Tkt, it decays exponentially: 



G v (x) 



=< e #(*)/2 e -i<«o)/2 



(6) 



The second factor in Eq.|5| is the gauge invariant 
Schwinger Green function of the spinon in a given con- 
figuration of a static gauge field: 

[GT v (x,t)] n =< ^ la (x,t)e l Io a ^ x)dx 4 lQ (0,0) > (7) 

The physical meaning of Schwinger gauge invariant 
Green function has been painfully sought. Remarkably, 
it is this function directly involved in the ARPES in the 
pseudogap regime in high temperature superconductor 

It is well known that the infrared divergence 
in the gauge fluctuation in Eq^ plagues the non- 
gauge invariant spinon Green function G s (x,t) =< 
ipiia(x,t)ipl la (0,0) >. As shown in it leads to 

divergent quasi-particle scattering rate I/77. However, 
the infra-red divergence has no physical meaning and 
was shown to be cancelled by the vertex correction in 



the gauge invariant transport time r fr Jig ]. Very simi- 
larly, this divergence will disappear in the gauge invari- 
ant G % ™ (x, t) which directly determine the ARPES data. 
The motion of non-relativistic particle in random mag- 
netic field has been studied in Ref. |2J|. It was found 
that the particle's propagation in random magnetic field 
ceases to be coherent and G™ v (x, t) develops branch cut 
singularities. Therefore quasi-particle picture completely 
breaks down which is a hallmark of non-Fermi liquid pic- 
ture. Studying massless Dirac fermion in random mag- 
netic problem is a much more complicated problem, we 
leave it to future publication p3). Instead, we can look 
at a much simpler problem of Dirac fermion interact- 
ing with dynamic fluctuating gauge field. This simpler 
problem has Lorenz invariance at T = which dictates 
Gl nv (k) ~ ik^^/k 2 - 2r i with k = (p,iu n ). It can be 
explicitly written as 
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n + Vfp x VAPy 
VAPy i^n - Vfp x 

(8) 



where Ei = (vfp x ) 2 + (vAPy) 2 and r\ is the anomalous 
dimension which was calculated to one loop in 1/N ex- 
pansion to be 77 - 0.27 in Ref. ||]. 

Taking analytic continuation to the real frequency and 
then taking imaginary part, we get the spectral function 
A(k, oj) = —ImG{k, ui + ie) of the electron near Ki at 
temperature T : 



A(k,uj) 



d 2 p 



sin(7r?7) 



VfPx 



(2tt)2 (£_£)a + £- 2(T) [lu 2 -E 2 }^-i 

x[6(uj - E P ) - 6{-uu - E P )] (9) 



with the similar expression at 
It can be shown that 



A(k,uj 
A(k,u) - 



0) 



1(7777) |t 



11 + 27, 



47r(fc 2 +£- 2 (T)) 



-1+217 



i(nr])\c 



(10) 



Precise energy distribution curve (EDC) ( at fixed k) and 
momentum distribution curve (MDC) ( at fixed lj ) for 
£(T) ~ 10 are obtained numerically and drawn in Fig. 2. 
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Fig 2: Energy distribution curve (EDC) ( at fixed k) and Mo- 
mentum distribution curve (MDC) ( at fixed w ). 

It is easy to see that just as random static gauge field, 
the dynamically fluctuating gauge potential also destroys 
the coherent quasi-particle picture. The spinon in the 
pseudo-gap regime loses its identity and replaced by col- 
lective excitations, in contrast to its counter-part in the 
superconducting state. The Volovik effect is short-ranged 
as shown in its effect is far less dramatic than that 
of the random gauge field which completely destroys the 
quasi-particle picture. Our interpretation of the ARPES 
data is in sharp contrast with that of Z2 gauge theory 
p5j. In this theory, the gauge field is absent due to the 
generation of double strength hc/e vortices which leads to 
spin charge separation in the long length scale at T = 0. 

From the long time behavior of G(r) ~ G" lv {t) ~ 
l/r 2+2? ', we find the local density of state (DOS) N{uj) - 
lo 1+2ti . However, as shown in |15[ the Volovik effect will 
generate finite DOS at zero energy N(ui = 0, T) ~ ^fnj ~ 
So the EDC curves in Fig. 2 should approach 
a small finite value at very small frequency due to the 
Volovik effect. 

A few caveats should be pointed out before comparing 
Fig.l with the available ARPES data in the node direc- 
tions: (1) The mean field approximate spin-charge sepa- 
ration Eq.|| is only valid at vortex hydro-dynamic regime 
( Fig.l ). It is only 30 — 40K above Tkt consistent with 
the recent experiments by Corson et al. Raising tem- 
perature higher will move into QC regime. There is no 
spin-charge separation at T — 0. A different calculation 
is needed to calculate ARPES spectrum in QC and QD 
regimes [Elj (2) Conceptually, the gauge invariant Green 
function should be calculated in a static gauge field in- 
stead of in a dynamic gauge field. However, we do not 
expect the technical calculation in static gauge field will 
lead to dramatically different EDC and MDC than those 
in Fig. 2, because unlike the gauge non-invariant Green 
function, the gauge invariant one is free of infra-red diver- 
gences. Therefore, we expect Fig. 2 qualitatively describe 
the electron spectral weight in the vortex hydro-dynamic 
regime if the phase fluctuation scenario is the correct one. 
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nous referees for careful reading and helpful comments 
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